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Thermal conductivity measurements were performed on single crystal samples of the supercon- 
ducting filled skutterudite compounds PrOs4Sbi2 and PrRu4Sbi2 both as a function of temperature 
and magnetic field applied perpendicular to the heat current. In zero magnetic field, the low 
temperature electronic thermal conductivity of PrRu4Sbi2 is vanishingly small, consistent with a 
fully-gapped Fermi surface. For PrOs4Sbi2, however, we find clear evidence for residual electronic 
conduction as the temperature tends to zero Kelvin which is consistent with the presence of nodes 
in the superconducting energy gap. The field dependence of the electronic conductivity for both 
compounds shows a rapid rise immediately above H c i and significant structure over the entire vortex 
state. In the fully gapped superconductor PrRu4Sbi2, this is interpreted in terms of multi-band ef- 
fects. In PrOs4Sbi2, we consider the Doppler shift of nodal quasiparticles at low fields and multiband 
effects at higher fields. 



Fundamental research into superconductivity is ex- 
tremely important because of the potentially revolu- 
tionary technology were this phenomena available for 
widespread commercial use. With improvements in ma- 
terials science, new materials that are superconducting 
continue to be discovered and experimental studies on 
them shed light on where our knowledge about super- 
conductivity can be enhanced. The filled skutterudite 
family of materials is one such example, exhibiting many 
features that suggest in some variants its superconductiv- 
ity may be conventional, whilst in others unconventional. 
A comparison between the superconducting properties of 
each type of material is certainly instructive. 

Some of the more unusual properties have been seen 
in the heavy- fermion superconductor PrOs4Sbi2. For ex- 
ample, in angle-resolved magneto-thermal conductivity 
experiments jl] , a change in the symmetry of the small 
anisotropy in the conductivity is interpreted as evidence 
for a multiphase superconducting phase-diagram. In this 
case, switching between a phase that has two nodes in 
the superconducting order parameter to one which has 
four nodes, as the magnetic field is increased. More re- 
cently, the thermal conductivity has been measured to 
even lower temperatures 0, Surprisingly, the mea- 
surements show an absence of any residual electronic con- 
duction in zero field, which is inconsistent with the earlier 
measurements in either superconducting phase. More- 
over, the large magnetic field dependence is interpreted 
as resulting from multi-band effects through analogy with 
well-known multiband superconductors such as MgB-2 [J] . 

On the other hand, the isostructural but non-heavy 
fermion compound, PrRu4Sbi2, appears to be a conven- 
tional s-wave superconductor with exponential temper- 
ature dependencies observed in specific heat [H, 0] and 
superfluid density Q. 



In this Letter, we present low temperature thermal 
conductivity measurements that provide evidence for a 
fully-gapped order parameter in PrRu4Sbi2 and a nodal 
superconducting order parameter in PrOs4Sbi2. We 
also demonstrate that the field dependence of the fully- 
gapped superconductor is consistent with multiband su- 
perconductivity. In the nodal superconductor, the field 
dependence immediately above H c \ is consistent with 
being due to the Doppler shift of nodal quasiparticles, 
with additional multiband effects occuring at fields above 
H C 2/2. Finally, we explore the possibility that multi- 
band superconductivity is a generic feature of filled- 
skutterudite superconductors. 

The single crystal samples were oriented using Laue x- 
ray backscattering. In the case of PrOs4Sbi2, the sample 
was then polished to a cuboid of dimension (600 x 100 x 
100) /im. For PrRu4Sbi2 the sample used was as-grown 
with dimensions (2500 x 188 x 380) /im. In both cases, 
four silver wires were then attached using high purity 
indium solder. 

The thermal conductivity k was measured using a sin- 
gle heater-two thermometer method. The heat current 
was supplied along the a-axis direction and the magnetic 
field applied perpendicular to this. The measurements 
were made in a dilution refrigerator by varying the tem- 
perature from 0.04 K to > 0.7 K at fixed magnetic field. 
The samples were field-cooled by cycling to T > 2 K be- 
fore changing the field. The error in the absolute value of 
the conductivity is estimated to be approximately 10%. 
The relative error between temperature sweeps at differ- 
ent fields is of order 1%. 

Zero Magnetic Field: In Fig. [T] the temperature depen- 
dence of the zero-field thermal conductivity is plotted for 
both materials. Since the measured quantity is the total 
thermal conductivity, it is necessary to separate the con- 
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FIG. 1: Thermal conductivity divided by temperature T ver- 
sus T 2 in zero field for PrOs4Sbi2 and PrRu4Sbi2. The lines 
are linear fits to the low temperature data (T < 0.15 K) which 
are then extrapolated to the T = axis. 



tributions from electrons and phonons. Using the well- 
established model based on kinetic theory (see, for ex- 
ample, 8]), in the low temperature limit, the electronic 
contribution is linear in temperature, while phonons are 
cubic in temperature. Fitting the measured conductivity 
(k) to the form: 



T T 



(1) 



and extrapolating the fit to zero temperature we obtain a 
value for the residual electronic conductivity divided by 
T, kq/T. The coefficient f3 in the above expression rep- 
resents the phonon contribution in this simple analysis. 
Fitting the data below T = 0.15, the values obtained for 
each material are 0.46±0.07 mW/K 2 cm and 0.058±0.007 
mW/K 2 cm for PrOs4Sbi 2 and PrRu4Sbi2 respectively. 

In PrRu4Sbi2, the value of kq/T is an order of mag- 
nitude smaller than that measured in PrOs4Sbi2- In in- 
stances where the extrapolated value is very small, the 
simplicity of the model used to separate electrons and 
phonons is exposed. Improved fitting accuracy has been 
demonstrated using an empirically determined variable 
power law for the phonon contribution [loj ]. In such 
cases the extrapolated value for the residual linear term 
is always found to be reduced and in this case we find 
k /T = 8± 17^W/K 2 cm and a phonon exponent of T 16 . 
Within error bars such a value can be considered as zero 
indicating the absence of any residual electronic conduc- 
tivity and therefore consistent with a fully-gapped super- 
conducting state. 

In contrast, the observation of a significant finite value 
for the zero-temperature extrapolation of the linear elec- 



tronic conductivity in PrOs4Sbi 2 is incontrovertible ev- 
idence for nodes on the Fermi surface of the supercon- 
ducting order parameter. 

As a consistency check, a calculation of the magni- 
tude of the residual linear term assuming a simple d- 
wave symmetry of the superconducting gap Q and weak 
coupling can be made for both systems. In PrOs4Sbi2, 
we estimate k /T = 0.3 mW/K 2 cm. However, we note 
that such a 4-fold symmetry would only be consistent 
with the high-field superconducting phase observed by 
angle-resolved magneto-thermal conductivity measure- 
ments For PrRu4Sbi2, this simple calculation would 
give an even larger value because of the smaller gap mag- 
nitude (assuming a similar Fermi velocity), which is con- 
sistent with very small extrapolated value being inter- 
preted as evidence for a fully gapped Fermi surface. 

In spite of the more than order of magnitude differ- 
ence in the electronic conductivity, the apparent phonon 
contribution in each material is very similar, as might be 
expected given the identical structure of these materials. 

The observation of a linear term in PrOs4Sbi2 in this 
study directly contradicts the observation of other low 
temperature thermal conductivity measurements where 
the data extrapolates to zero with a rapid temperature 
dependence that exceeds T 3 0, Q • Rapid temperature 
dependencies such as these are known to occur when 
electron-phonon decoupling hampers the measurement of 
the intrinsic electronic conductivity 11| and is thought 
to occur in other superconducting systems [H, El- In 
CeCoIn 5 in particular, early measurements indicated a 
large T 3A temperature dependence and zero electronic 
residual conductivity [ljj]. More recent measurements 
show a large residual electronic conductivity and a lower 
T 2 temperature dependence [l5j |. In the study reported 
here, the use of an identical technique to measure each 
material with identical contacts and yet dramatically dif- 
ferent but nevertheless consistent results is compelling. 
Magnetic Field Dependence: The temperature depen- 
dence of the thermal conductivity for different magnetic 
fields applied perpendicular to the direction of the heat 
current is shown in Fig. [5] for PrOs4Sbi2 and Fig. [5] for 
PrRu4Sbi2- 

The data for both materials show features that are 
qualitatively similar, although there are some distinct dif- 
ferences. In zero field, the data have the steepest temper- 
ature dependence, which in PrOs4Sbi 2 rises to a plateau 
above T = 0.6 K. Measurements at higher temperature 
show this to be a peak [X71 | . In both materials, as soon 
as a magnetic field is applied the T-dependence is sup- 
pressed by a dramatic increase in k/T at the lowest tem- 
peratures. For subsequent higher fields, the curves at 
low temperatures then remain approximately parallel as 
the field increases by orders of magnitude and the sys- 
tem enters its field-induced normal state. We note in 
passing that the data for PrOs4Sbi2 at low fields show 
a striking similarity to the behaviour observed in ultra- 



3 



1 1 1 1 1 1 1 






— 


* * 

- 


8 


- 


• • • . 












Aon ^On A n a 


- j^"** * .' PrOs 4 Sbi2 


B ^.O 1 


• 

• 


♦ LOT 


• 


• 0.4 T 


• 

• 

• 


= 0.2 T 


• 

• 


- 0.1 T 


o 55 mT 


i i i i i i 


• zero field 



10 



0.0 



0.5 

Temperature [K] 



1.0 



FIG. 2: Thermal conductivity divided by temperature T 
versus T for PrOs4Sbi2 at different magnetic fields. The heat 
current is applied along the a-axis of the single crystal sample. 
The magnetic field is applied perpendicular to this direction. 



E 
o 

CM 

E, 

t- 

"2 



0>— 
0.0 



1^ 



0.5 

Temperature [K] 



PrRu4Sbi2 

» 2.0 T 

• 0.15 T 
' 0.10 T 
° 30 mT 
' 7 mT 

• zero field 



1.0 



FIG. 3: Thermal conductivity divided by temperature T 
versus T for PrRu4Sbi2 at different magnetic fields. The heat 
current is applied along the a-axis of the single crystal sample. 
The magnetic field is applied perpendicular to this direction. 



pure YBa 2 Cu 3 07 UM- 

In order to explore the magnetic field dependence more 
clearly, the value of the temperature dependence extrap- 
olated to T — K at each field is extracted, normalised 
to the normal state value, and then plotted in Fig. [4] as a 
function of magnetic field normalised to the upper criti- 
cal field {H C 2). Values used for the upper critical field are 
H c2 = 2.0 T and 0.2 T for PrOs 4 Sbi 2 and PrRu 4 Sbi 2 , re- 
spectively. In the simplest fully gapped superconductor, 
the field dependence of the residual electronic conduction 
rises exponentially with magnetic field on a scale se t by 
H C 2, as has been observed experimentally in VaSi [19j |. 
In contrast, Fig. Q] shows the field dependence for both 
materials rises rapidly at low fields. On this reduced mag- 
netic field scale, the increase is considerably more rapid 
in PrOs4Sbi2 than in PrRu4Sbi2. Futhermore, over the 
entire field range both materials show appreciable struc- 
ture. Such behaviour is reminiscent of that seen in multi- 
band superconductors such as MgB2 Q and NbSe2 [l9j 
and has been suggested already for PrOs 4 Sbi2 0,. Since 
the zero field behaviour has already established the fully 
gapped nature of the superconducting order parameter 
in PrRu4Sbi2, this is the most likely explanation in this 
material. 

In PrOs4Sbi2, the picture is less clear because of the 
presence of nodal quasiparticles. The principal effect of 
magnetic field on superconductors with a nodal order 
parameter is the Doppler shift of the nodal quasiparti- 
cle energy spectrum through coupling to the superfluid 
flow around magnetic vortices. This leads to an increase 
in the density of states and associated properties. The 
original semi-classical treatment by Volovik [l8| has been 



extended to explain the effect on thermal transport by 
Hirschfeld and Kubert with the increase dependent 
on the normal state scattering rate and given by 



/c(0, H) kq 
T ~ Y 



?yl + p 2 — sinh 1 p 



(2) 



The parameter p = ^/3<I>o7 /ahvp\/H is related to the 
impurity bandwidth (7) and the normal state scattering 
rate (r), 7 ~ T 2 . In the inset to Fig 31 we plot the nor- 
malised increase in thermal conductivity at T — 100 mK 
with magnetic field for both the data presented here and 
for the earlier measurements Q. In each case, the initial 
increase is fitted to extract values for p. The values ob- 
tained are p — 90 for this study and p — 183 for the data 
from 0] • We note that there is some variability in exactly 
what range of data to fit since at higher fields there is a 
crossover to a field independent region. We have chosen 
to fit the data up to 80% of the plateau value conduc- 
tivity. Adjusting the range in each case causes a shift in 
the values of p by approximately 10%. Nonetheless, the 
important point is that the values for p in each case differ 
by a factor of two, which implies that the normal state 
scattering rates and, hence, normal state conductivity 
should vary by a factor of four. Comparing our normal 
state conductivity value, k (H — 2T)/T = 8 mW/K 2 cm, 
with that from RefQ], k (H = 2T)/T = 2 mW/K 2 cm, 
we see that quantitatively this is exactly the case. The 
rapid increase with magnetic field therefore scales with 
the normal state scattering rate in a manner that is en- 
tirely predicted by a Doppler shift of the quasiparticle 
energy spectrum. Unfortunately, and as has been es- 
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tablished already from similar measurements in cuprates 
the absolute value of the normal state scattering 
rate obtained from these fits is not in agreement with 
that from the normal state transport properties. This 
most likely reflects the limitations of the current theory. 

At higher fields in PrOs4Sbi2, the conductivity ex- 
hibits a plateau and then increases in a super- linear fash- 
ion up to the normal state value at H c i- This aspect 
is again reminiscent of the behaviour seen in multiband 
superconductors. Consequently, we postulate that this 
high-field dependence results from quasiparticles associ- 
ated with a second, larger superconducting gap covering 
a separate Fermi surface sheet. The resulting picture for 
PrOs4Sbi2 is that it has one band with a nodal order pa- 
rameter and another distinct band that is fully gapped. 
This scenario is appealing for a number of reasons. First, 
if we postulate that the gap maximum on the nodal band 
is somewhat smaller than the gap maximum on the fully 
gapped band, then our estimate for the value of the resid- 
ual electronic conductivity would be revised upwards and 
closer to the value we measure. Second, the rapid rise of 
the conductivity at low fields would now be considered 
on a scale set by a presumably lower critical field for the 
band with a nodal gap. This may provide a more reason- 
able absolute value for the normal state scattering rate. 
Third, it may provide a consistent picture of the angle re- 
solved magneto-thermal conductivity measurements [H. 
If we postulate a degree of anisotropy to the fully gapped 
system with a minimum orthogonal to the nodes of the 
other gap, then high field or high temperature measure- 
ments would pick up a four-fold symmetry due to the 
nodes on the gap of one band and the anisotropy on the 
fully-gapped other band. At low fields and temperatures, 
only the anisotropy due to the nodal sheet would survive 
with quasiparticle excitations on the (anisotropic) fully- 
gapped band exponentially suppressed. Such a system 
would thus have a superconducting state that is not only 
multi-band, but also multi-symmetric. The reason be- 
hind why one order parameter should have nodes and the 
other fully-gapped may be related to the heavy-fermion 
character of particular sheets of the Fermi surface. 

In conclusion, we have measured the low temper- 
ature thermal conductivity of the filled skutteruditcs 
PrRu4Sbi2 and PrOs4Sbi2 as a function of transverse 
magnetic field. Measurements in zero field indicate that 
in PrRu4Sbi2, there is no evidence of any residual con- 
duction, consistent with a fully gapped superconducting 
state. Conversely, the extrapolated data in PrOs4Sbi2 
shows a distinct linear contribution consistent with low 
energy electronic excitation and therefore with a super- 
conducting gap with nodes on the Fermi surface. A rapid 
increase of the conductivity with magnetic field is ob- 
served in both materials. In PrRu4Sbi2, the increase 
in conductivity is less rapid but still above what one 
would expect for a simple isotropic gapped superconduc- 
tor. This increase is attributed to multiband supercon- 
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FIG. 4: Magnetic field dependence of the extrapolated 
zero temperature thermal conductivity of PrOs4Sbi2 and 
PrRu4Sbi2. The lines are guides to the eye in each case. 
The inset shows the field dependence of the conductivity nor- 
malised to the zero-field valued for PrOs4Sbi2 both from this 
study and from earlier work 2j. The lines result from fitting 
to a semi-classical theory [20] based on a Doppler shift of the 
nodal quasiparticle spectrum. 

ductivity. In PrOs4Sbi2, the increase is consistent with 
nodal quasiparticle excitations coupling to the superfluid 
flow around magnetic vortices. Moreover the magnitude 
of this increase is found to scale appropriately with the 
normal state scattering rate when compared to existing 
data for measurements on samples with shorter quasipar- 
ticle lifetimes. At higher fields, the superlinear increase 
in conductivity is attributed to quasiparticle excitations 
associated with a large fully-gapped band on a separate 
sheet of the Fermi-surface. If such a scenario is correct, 
it is possible that multiband superconductivity may be 
a generic feature of superconducting filled skutterudite 
systems. As this study demonstrates, this appears to be 
the case for both PrR.U4Sbi 2 and PrOs4Sbi 2 . 
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